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ABSTRACT 

The outer disk of the Large Magellanic Cloud (LMC) is studied in order to 
unveil clues about its formation and evolution. Complementing our previous 
studies in innermost fields (3 kpc<R<7 kpc), we obtained deep color magnitude 
diagrams in 6 fields with galactocentric distances from 5.2 kpc to 9.2 kpc and 
different azimuths. The comparison with isochrones shows that while the oldest 
population is approximately coeval in all fields, the age of the youngest popula- 
tions increases with increasing radius. This agrees with the results obtained in 
the innermost fields. Low-resolution spectroscopy in the infrared Ca II triplet 
region has been obtained for about 150 stars near the tip red giant branch in 
the same fields. Radial velocities and stellar metallicities have been obtained 
from these spectra. The metallicity distribution of each field has been analyzed 
together with those previously studied. The metal content of the most metal- 
poor objects, which are also the oldest according to the derived age-metallicity 
relationships, is similar in all fields independently of the galactocentric distance. 
However, while the metallicity of the most metal-rich objects measured, which 
are the youngest ones, remains constant in the inner 6 kpc, it decreases with 
increasing radius from there off. The same is true for the mean metallicity. Ac- 
cording to the derived age-metallicity relationships, which are consistent with 
being the same in all fields, this result may be interpreted as an outside-in for- 
mation scheme in opposition with the inside-out scenario predicted by ACDM 
cosmology for a galaxy like the LMC. The analysis of the radial velocities of our 
sample of giants shows that they follow a rotational cold disk kinematics. The 
velocity dispersion increases as metallicity decreases indicating that the most 
metal-poor/oldest objects are distributed in a thicker disk than the most metal- 
rich/youngest ones in agreement with the findings in other disks such as that of 
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the Milky Way. They do not seem to be part of a hot halo, if one exists in the 
LMC. 

Subject headings: galaxies: abundances - galaxies: evolution - galaxies: kinematics - 
Magellanic Clouds - galaxies: stellar content 
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Introduction 



The formation and evolution of stellar disks in galaxies are still not well understood 
proc esses in spite of the observation al and theoretical effort developed in the last years 



see 



van der Kruit fc Freeman 



2011 



for a review). Simulations based on ACDM predict 



that the inner parts of disks form first and then they grow up to reach their present-day 



size (e.g. 



White fc Frenk 



1991 



Mo et al. 



1998 



Roskar et al. 



2008|). According to this 



scenario, so-called inside-out, it is expected that the youngest populations are located in 
the outermost regions. Therefore, to understand the formation and evolution of disks, it 
is important to have information about the spatial distribution of their different stellar 
populations in a large range of galactocentric distances. In particular, their outskirts are 
key since they might retain a fossil record of the disk formation due to their long dynamical 
timescales. 

From an observational point of view, the study of the outer disks is difficult due to their 
low-surface brightness at the level of, or even bellow, the background. The behaviour of the 
star formation r ate with radius can be studied w ith surface photometry in galaxies within 



a few Mpc (e.g. 



de Jong 



1996 



Taylor et al. 



20051 ). In general, t hese investigations found 



that massive early-type spirals are bluer in their outer disk (e.g. 



Trujillo fc Pohlen 



20051 ) 



in agreement with the inside -out scenario although these d isks are mainly populated by 



old (~10 Gyr) old stars (e.g. 



Sanchez-Blazquez et al. 



201 lh . Disks i n late-type spirals and 



de Jong 



1996 



Gadotti fc dos Anjos 



irregular galaxies have more puzzling behaviors (e.g. 

1 Based on observations collected at the European Southern Observatory, Chile, within 
the observing programs 074.B-0474 and 082.B-0900 

2 The National Radio Astronomy Observatory is a facility of the National Science Funda- 
tion operated under cooperative agreement by Associated Universties, Inc. 

3 Adjoint Professor 



20011 ; iTaylor et al.l 120051 ; iHunter &: Elmegreenl 120061 ). Some of them follow the same trend 
observed in early-type spirals. However, others become redder outward which is interpreted 
as the youngest populations being located in the innermost regions, since no galaxy 



formation model predicts 
in the way pr e dicte d by 



positive meta l 



Stinson et al. 



icity gradients. This suggests an outside-in scenario 



( 120091 ) for less massive galaxies and observed by 



Hidalgo et al.l ( 120111 ) for less massive (M<10 9 M ) galaxies. 



Nearby galaxies, in which their stars can be resolved, offer the opportunity to study 
galactic disks in more detail. The observations in the two most massive spirals in the Local 
Group, M31 and the Milky Way, agree quite well with the predictions of the inside-out 



scenario. M31 does not show a clear metallicity gradient (jlbata et al. 



20051 ) while in 



the Milky Way a smoot 



i gradient is observed which was more sharp in the past (e.g. 
201 ll ). There is another spiral in the Local Group, M33, of later type 
and much less massive than the previous ones. This galaxy shows a truncation in its 



Carrera fc Pancino 



radial bright ness surface profile 



galaxies (e.g. 



Trujillo fc Pohlen 



( Ferguson et al. 



20071 ) which is also observed in other 



20051 ). Although still not well un derstood, this trunc ation 



can be originated during the collapse of the protogalactic cloud (Ivan der Kruit 



19871 ) or 



can be related to different star formation thresholds at different galactocentric distances 



Schave 



( j Williams et al. 



2004 ). T he scale-length of M33 increases with time with in the break radiu s 



2009) and steepens and even reverse from there out (IBarker et al 



20071 ). 



These authors interpreted this result in the context of the inside-out scenario: the inner disk 
would have been grown from the inner regions, while secular processes radially redistribute 
stars causing an inflection point in the stellar age gradient at the break radius. 

There is a fourth disk galaxy in the Local Group, the Large Magellanic Cloud (LMC). 
Sometimes considered an irregular, it is more properly classified as a Barred Magellanic 
Spiral (SBm). This kind of objects are frequent through the Universe and are characterised 
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by a mass between that of large and dwarf galaxies, and the presence of a bar and a disk 
containing a single spiral arm. The proximity of the LMC allows us to investigate it using 
both spectroscopy and photometry of individual stars. Therefore, we can obtain information 
about the age and/or chemical abundances spatial distribution and about kinematics. 
Moreover, the LMC contains the only resolved bar in the Local Group and it forms an 
interacting system together with the Small Magallanic Cloud and the Milky Way, a l thoug h 
their history of past interactions are still not well understood (e.g. 



Bcsla ct al. 



2007 



2010|). 



2000 



2008 



Castro et al. 



e^ 



The LMC inner disk^ 
stellar population 



(R < 4.5 



Bertelli et al 



200 



Harris fc Zaritsky 



(jMajewski et al. 



2009 



ipc) has been the subject of most work about its 



1992 



O 



Smecker-Hane et al. 



sen 



1999 



2002 



Holtzman et al. 



Cole et al. 



2005 



1999 



Cole et al 



Pompeia et al. 



2009). H oweve r, the LMC disk extends at least to 12-14 kpc 



Saha et al. 



2010). The lack of works in the outer LMC do not allow 



to study radial trends in deta il. Studies covering large galactocentric distances use stellar 



clust ers as test particles (e.g. 



Olszewski et al 



1991 



Grocholski et al 



2006 



Sharma et al 



20101 ). They did not find neither age nor metallicity gradients. The spatial distribution of 
asymptotic giant branch stars (AGB) in the field of the LMC has been also studied using the 
ratio between the number of carbon and oxygen-rich stars, traditionally denoted as C/M. 
Two works have used this indicator reaching contradicting conclusions for the inner disk. 



Cionil ( 120091 ) found that the C/M ratio increases rapidly with increasing distance which is 



interpreted as a met allicity gradient wi th the outskirts being more metal-poor. However, 



using the same data 



Feast et al. 



(120101 ) found that the C/M ratio is nearly constant in 



the inner 4 kpc, from there off the ratio decreases. Since carbon-rich stars are rare in 



4 In this paper we call "inner disk" the region within R ~ 4.5 kpc, and "outer disk" the 
region from that distance off. Note t hat at the distance of the LMC, 1 kpc corresponds to 
1?12 on the sky (Ivan der MarellboOli ) 
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old populations this result is interpreted as the outer disk stars are older. Unfortunately, 
high-resolution spectroscopy studies at different galactocentric distances, which would allow 
investigation of the possible presence of abundance gradients of some key chemical species, 
are not available. 

Our group is involved in a major effort to stud y the LMC using both photometric 



and spectroscopic observations. 



Gallart et al. 



(12008 



hereafter 



; 



aper I) studied four fields 



located at different galactocentric distances (3 kpc <R< 7 kpcjj. They found that the age 
of the youngest populations in each field increases with galac tocentric distance. In fact, a 



more detailed analysis of these fields by 



Meschin et al. 



(201 if) found a gr adual decrease of 



Carrera et al. 



(j2008al hereafter Paper 



the ratio of young to old populations with radius. 
II) derived ages and metallicities for about 100 red giant branch (RGB) stars confirmed as 
LMC members from their radial velocities in each of these fields using the infrared Ca II 
triplet (CaT) lines. The three inner fields have a similar mean metallicity, but that of 
the outermost one is a factor of two more metal-poor (see Paper II for details). This is 
explained by the lack of the youngest stars in the latter field which are also more metal-rich. 
These results point toward an outside-in formation of the LMC disk, on contrary to the 
inside-out scenario usually predicted for a galaxy with a mass similar to tha t of the LMC. 



The observed behaviour i s more like that expected for less massive galaxies ( iStinson et al 



2009 



Hidalgo et al. 



201lh . 



In this paper, we extend our study of the LMC disk at larger galactocentric distances 
(R>7 kpc). For this purpose we have obtained deep CMDs and low resolution spectroscopy 



5 The galactocentric distances used in Paper I are different from those used here and in 
Paper II, because they have been calculated assuming different centres. While the dynamical 
center of the Carbon stars was used here, which better represents the center of the older 
populations, paper I used the center of the Bar. The distance between both is about 1 kpc. 
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for 6 new fields spread over the LMC disk. Their spatial distribution and color-magnitude 
diagrams (CMD) are described in Sections [2] and [21 respectively. The spectroscopic 
observations and data reduction, together with the tools used to derive radial velocities, 
metallicities, and ages are described in Section HI The metallicity distribution and age- 
metallicity relation (AMR) of each field are presented in Sections 17.11 and 17.21 respectively. 
The kinematics of our RGB sample is analysed in Section 17.31 Finally, the implications of 
the main results of this work are discussed in Section [8] and summarised in Section [91 



2. Spatial Distribution of Observed Fields 



We have selected 6 fields in addition to the ones presented in P apers I and II 



van der Marell (120011 ) 



'hey 



are shown in Figure [T] together with the near infrared isopleths from 
for reference. Two of the new fields presented in this paper are situated at either side of 
the LMC bar, in the East- West direction, toward the SMC (West) and the Galactic center 
(East). With them we aim to investigate the possible influence of the SMC in the LMC's 
disk evolution. The remaining two fields are located northward of the bar, as those studied 
in Paper II, but further out than them. These fields were selected to study the outermost 
regions of the LMC disk and are key to investigate the presence of a metallicity gradient. 

Equatorial -right ascension and declination- coordinates over the celestial sphere for 
observed fields, named according to their right ascension and declination as in Paper I, are 



listed in Table [U Angular -dista nce and position ang 



using the equations described by Ivan der Marell ( 120011 ) consideri ng the dynamical center o f 



e- coordinates have been computed 



van der Marel et al 



20021) 



the carbon stars as origin (a^ooo = 5 h 2T^6 and #2000 = — 69?87; 
are also listed in Table [TJ Distances computed in this way might not be actual physical 
distances because they do not take into account the viewing angles of the disk: inclination 



i, and the major-axis position angle, P.A. ma ~, which are needed to deproject the observed 
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morphology. We deproject the angular distances of o ur fields iL.ni following 



assuming i = 34? 4 ± 6?2 and P. A. 



rnaj 



189?3 ± 1?4 ( Ivan der Marel fc Cioni 



Cionil ((20091), 



20011 ). These 



values, listed in the last column of Table [Q will be used throughout this work. 



3. Color-Magnitude Diagrams 

V- and /-band images for the northern fields, and B- and i?-band images for the 
others, were obtained with WFI at the 2.2 m telescope at La Silla Observatory, Chile. The 
WFI camera is a mosaic of 4 x 2 of 2048 x 4096 CCDs covering a total area of 34' x 
33' with a scale of 0.238 pixel -1 . The images were bias and flat-field corrected in the usual 
way with the IRAFlj MSCRED package. Profile- fitting photometry was obtained with the 



DAOPHOT /ALLFRAME suite of c odes ( 



Stetson 



system using observations of several IStetsonl ( 120001 ) fields obtained in the same run. 



19941 ) and calibrated to the standard 



\{B — J)o,Mj] CMDs of the western (towards the SMC) and eastern (in the opposite 
direction) fields are shown in the bottom and top panels of Figure |2j respectively. [{V — J) 
Mi] CMDs of the northern fields are plotted in Figure [3j Although we postpone to a 
forthcoming paper the detailed analysis of these CMDs, a set of isochrones has been 
overimposed on them with the purpose of intercomparing fields. We us ed the same set o 



overs hooting; isochrones as in Paper I selected from the BaSTI library (jPietrinferni et al. 



20041 ) . Their metallicities follow approximately the AMRs derived in Paper II (see section 



a. 



All fields studied, even those at a distance of ~9 kpc, have a well populated old 



6 IRAF is distributed by the National Optical Astronomy Observatories,which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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main-sequence turn-off region. This, in agreement with the results obtained in Paper I, 
suggests that star formation started at about the same time in all of them. An alternative 
is that these stars were formed in the inner disk and then they have migrated to the 
outskirts. All fields contain an important amount of intermediate-age LMC populations. 
The bulk of stars in the innermost field studied, LMC0422-7100, have ages older than 2 
Gyr while most of the stars in the remaining fields have ages older than 5 Gyr. In all fields 
a number of stars overlapping the 1 Gyr can be observed, which number increases as the 
global nu mber of the stars in the field increases. The presence of this population was also 



noted by 



Gallart et al. 



(120041 ) in a field located at ~7 kpc and interpreted as a burst of 
star formation 1 Gyr ago. However, since they are present in all studied fields and their 
number is correlated with the density of stars in each field they may be blue stragglers. As 
in the case of Paper I the comparison with isochrones show an increase of the age of the 
intermediate-age populations present in each field with galactocentric distance. There are 
not clear differences with azimuth since CMDs of fields at similar distances but different 
azimuths have a similar morphology. 



4. Spectroscopic Observations, Data Reduction and Radial Velocities 

As in Paper II, the stars observed spectroscopically were selected on the upper part of 
the RGB trying not to exclude metal-poor objects on the blue ridge but avoiding inclusion 
of red-clump and red supergiant objects. Stars below the tip of the RGB were given the 
highest priority of observation and brighter objects were only targeted when no others 
were available. The observations were carried out in service mode on ESO period 82 with 
FLAMES mounted at VLT UT2. We used the low-resolution grating LR8 which covers the 
wavelength range between ~ 8200 to ~ 9400 A with a resolution of R ~ 6500, in order to 
accurately measure equivalent widths of the CaT lines as well as radial velocities. 
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The ESQ pro vides reduced data obtained with the GIRAFFE data-reduction pipeline 



flMelo et al. 



20091 ). This pipeline performs the bias, flatfield and scattered light corrections, 
find, trace and extract the spectra and obtain the wavelength calibration based on daytime 
calibration exposures. The only reduction step that this pipeline does not perform is the 
sky subtraction. To do this, we followed a similar p rocedure as in Paper II but with some 



modifications as described by 



Battaglia et al. 



( 120081 ). In brief, for each exposure we average 



the spectra of all fibers placed on the sky to obtain a high signal-to- noise ratio (S/N) sky 
spectrum. This sky spectrum, and the objects spectra are separated into two components: 
continuum and line spectra. To obtain the continuum we fitted a polynomial rejecting the 
contribution of both sky and object lines. The line spectrum is obtained by subtracting 
the continuum. The sky-line component is cross correlated with each object-line to put 
both in the same wavelength scale. This also provides an additional test of the wavelength 
calibration. After updating the wavelength calibration, the sky- and object-line components 
are compared looking for the scale factor that minimises the sky line residuals over the 
whole spectral region. The object-continuum component is wavelength updated and added 
back to the sky subtracted object-line spectrum. Finally, the sky-continuum is subtracted 
assuming that the scale factor is the same than for the sky-line component. The result is 
that the residuals of the sky lines are lower than ~3%. 

Since 2 or 3 exposures were obtained for each configuration we combined all the spectra 
of the same object, thus improving the S/N and minimising the contribution of cosmic rays 
and bad pixels. Finally, each object spectrum is normalised. 

Radial velocities are key to reject foreground objects, but they also provide valuable 
information about the LMC kinematics (see Section [773]) . However, radial velocity standard 
stars were not observed, so we used one of our targets with a very high S/N as template. We 
obtained the radial velocity of this star using the laboratory position of the three CaT lines 
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and other lines easily identified in its spectrum using the rvidlines task in IRAF. Then, the 
radial velocities of all stars were calculated by cross-correlation using the IRAF fxcor task, 
with this high S/N star as template. As in paper I, we conside red as LMC mem bers those 



stars with radial velocities in the range 170<V r <380 km s 1 (IZhao et al. 



20031 ). In total, 



we have determined ages, and metallicities for 147 stars confirmed as LMC members from 
their radial velocities in the 6 fields observed. Their magnitudes, CaT equivalent widths 
(see below), and radial velocities are listed in Table [2j 



5. CaT Equivalent Widths and Metallicity Determination 



The metallicity of each star was determined from the 



lines using the procedure described by 



Carrera et al 



equivalent widths of the CaT 



(120071 . hereafter Paper III). The line 



profiles were obtained by fitting a Gaussian plus a Lorentzian, which provides the best fit 
to line core and wings (see Paper III for details). Although the spectra had already been 
normalised, the position of the continuum was recalculated by a linear fitting to the mean 
values of each continuum bandpass defined in Paper III. The CaT index, denoted as ECa, 
is defined as the sum of the equivalent widths of the three CaT lines. 

The strength of the CaT lines depends on the temperature, gravity and chemical 
abundance of the stellar atmosphere. Since we are only interested on the abundance 
dependence, the temperature and gravity contributions must be removed. The procedure 
followed is described in depth in Paper III and we summarize here the most important 
details. For a given chemical abundance, the stars define a sequence in the Luminosity- E 
Ca plane when the temperature and/or gravity change. Although this sequence is not linear 
for the whole magnitude range covered by the RGB it is very close to linear when only stars 
in the upper two magnitudes of the RGB are considered as in our case. The non-linearity 
of the sequences in the Luminosity-S Ca plane is in general important for [Fe/H]<-3 
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( IStarkenburg et al.ll2010l ). However, the expected number of stars in that metallicity range 
is negligible based on the results described in Paper I. We safely assume therefore that the 
sequence defined by stars of the same metallicity in the Luminosity- E Ca plane is linear. 
We used Mj as luminosity indicator, because the CaT metallicity calibration based on Mj is 
less sensitive to age. In the M/-E Ca plane the sequences can be parametrised in the form 
ECa = W'j + (3jMj. W/ is called the reduced equivalent width. The slope is independent 
of the metallicity, whereas W/ changes with the chemical abundance. 

In Paper III, we obtained relationships between W j and metallicity in different 



metallicity scales, including the 



Carretta fc Grattonl (119971 hereafter CG97). The stellar 



metallicities o 



Carretta et al. 



;he o bjects studied in Paper II are calculated in this scale. Recently 



(120091 . hereafter C09) have updated their metallicity scale based on more 
homogeneous and higher resolution spectra of a larger sample of bright giants in 24 
Galactic globular clusters. They provide the relation between their new and old scales. The 
existence of this new metallicity scale for globular clusters as well as new determinations 
of metallicities for some open clusters used in Paper III from high-resolution spectroscopy 
(R~20000) motivate us to obtain a new relation between the strength of the CaT lines 
and metallicity. To do this we used the same cluster sample as in Paper III and listed in 
Table [31 The reduced equivalent width of each cluster in V and / shown in columns 2 and 
3 are the same of those used in Paper III. The metallicities used in Paper III, the new 
values calculated here and the corresponding references are listed in column 4, 5 and 6, 
respectively. 

The behaviour of W'y and W'j as a function of metallicity are shown in Figure HI We 
have used 25 clusters for the calibration in V and 23 for that in /. In both linear 
correlation is obtained with best fits given by: 
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[Fe/H] v = -3.62(±0.02) + 0.43(±0.01)W£, a v = 0.03 (1) 
[Fe/H] 1 = -3.38(±0.02) + 0.44(±0.01)W£ oj = 0.03 (2) 

In Paper III we determined that the linear sequence in the M/-E Ca plane is given by 
SCa = Wj — 0.611 x Mj. Therefore, using Equation [2] the metallicities of the LMC stars 
are given by: 



[Fe/H] = -3.38 + 0.44SCa + 0.27M/ 



(3) 



To obtain the absolute magnitudes of the LMC stars we assumed a dist ance modulus 



of (m — M)o = 18.5 ± 0.1 and the reddenings listed in Table [T] obtained from 
(119981 ). 



Schlegel et al. 



Since different metallicity scales are available in the literature, it is necessary to 
determine to what extent the metallicity distributions obtained here depend on the scale 
used. To do this we have calculated the metallicities of all stars observed in the northern 
fields (~ 400) in three different scales, which we display as follows. The left panel of 
Figure [5] shows the difference of metallicities obtained in Paper II and those obtained here 
as a function of the last ones. The new values are systematically more metal-poor than 
those derived in Paper II, and the difference increases towards lower metallicities. This is 



explained by the difference between the two globu 



case (for a detailed comparison see 



Carretta et al- 



ar clu ster metallicity scales used in each 



2009l ). It is interesting to note that all 
the CaT surveys performed until now in dwarf spheroidals companions of the Milky Way 



do no t find objects as metal-poor as those observed in the Galactic hal o (e.g. 



2006). The new scale seems to solve in part this discrepancy. However, 



Helmi et al 



Starkenburg et al 



(120101 ) demonstrated that the lack of metal-poor stars in CaT surveys is due to the way in 
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which the temperature and gravity contributions are taken into account. They proposed 
that extremely metal-poor stars are not detected because of the wrong assumption that 
stars of the same metallicity follow a linear sequence in the M/-E Ca plane (see also Paper 
III). On t he right panel of Figure [5] we have compared the difference of metallicities derived 



using the 



Starkenburg et al 



(120101 ) recipe to those obtained he re with our methq c 



before, our metallicities are systematically lower than those by IStarkenburg et al. 



As 



201(1 ) 



at low metallicities. However in this case the dispersion is much larger. We conclude that 
the differences between the three scales are only noticeable (A[Fe/H]>0.2 dex) for extreme 
values of the metallicity but that in the metallicity range which contains most of the LMC 
stars (i.e., -1.6<[Fe/H]<+0.3) they produce values within ±0.2 dex. 

The metallicity distributions on each field are shown in Figures |6] and [7J and will be 
discussed in extenso in Section 17.11 



6. Estimation of Stellar Ages 



The age-metallicity degeneracy in the position of the RGB in the CMD can be broken 
when either the metallicity or the age can be obtained independently of each other. Because 
of this, combining the metallicities of stars derived from spectroscopy with their position 
in the CMD allows us to estimate their age. In Paper II, a polynomial relationship was 
computed to derive stellar ages from their metallicities and positions in the CMD. F or 



this purpose, a synthetic CMD was created with IAC-star ( Aparicio fc Gallart 



200J), 



adopting BaSTI stellar evolution models with overshooting (jPietrinferni et al.ll2004j). The 



Carrera et al 



aforeme ntioned relation was obtained in Paper II for (V — I) and My while 
( ]2008bl ) derived a similar relation from the same synthetic CMD for (B — I) and Mj. We 
refer the reader to these papers for a detailed discussi on of the procedure used, uncertainties 
and tests of confidence levels 



Carrera et al 



( I2008bl ) confirmed that both relationships 
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produce similar results within the uncertainties. We have used the corresponding coefficients 
-those derived for (V — I) in the case of the North fields, and those for (B — /) for the 
remaining fields- obtained in these papers to derive ages for each program star. 

7. Analysis 

7.1. Metallicity Distributions 

The metallicity distributions for each field are shown in Figures [6] and [7J Because 
we used here a different metallicity scale than in Paper II, we have also recalculated the 
metallicity distributions of the fields analysed there (Figure [6]). The differences between 
the distributions described in Paper II and those obtained with the new scale are only 
appreciable for low metallicity stars, which are poorly sampled. The small number of 
objects (< 50) confirmed as LMC members from their radial velocity in the new fields 
studied, except for field LMC0422-7100 (see Table H]), do not allow us to perform a detailed 
analysis of these distributions. However, we can use their mean metallicity to investigate 
possible differences among them as well as metallicity gradients. In Paper II means and 
dispersions were obtained by Gaussian fitting to each distribution. This is correct for 
samples with a number of objects large enough (>50). But Gaussian fitting is not adequate 
for poorly sampled distributions. To overcome this problem we have adopted the following 
procedure. Half of the stars in each sample were randomly selected and the mean, standard 
deviation and 10th and 90th percentile of their metallicity distribution were computed. This 
procedure was repeated 10 4 times using random subsets each time. The values obtained in 
all of them for each statistical parameter have Gaussian distributions centered in the true 
mean, dispersion or percentile value. The width of each distribution provides an estimation 
of its uncertainty. The values thus obtained are listed in Table HJ where the fields are 
ordered by their distance to the center (see Section [2]) given in column 2. Uncertainties are 



-17- 



smaller for better sampled fields, as expected. The RGB sample studied in the LMC bar by 



Cole et al 



( 120051 ) has been also analysed in the same way and the results are also listed in 
Table HI Since their CaT index is equivalent to ours, as was demonstrated in Paper III, we 
used their index to calculate the metallicities from Equation |3j 

The mean metallicity values in the inner ~7 kpc are similar within the uncertainties 
but are lower than ~-0.8 dex beyond this radius suggesting the presence of a gradient. The 
computed 10th and 90th percentile values provide further insight into the possible presence 
of a gradient. They are listed in columns 5 and 6 of Table HJ and have been plotted in 
Figure E] as a function of radius. The 10th percentile value (filled symbols), account for 
the most metal-poor regime and have similar values in all disk fields and in the bar within 
the uncertainties. The 90th percentile values (open symbols) which account for the most 
metal-rich populations, are however, constant within the inner ~7 kpc, but decrease from 
this point on, indicating a metallicity gradient in the most metal-rich populations outwards 
of ~7 kpc. The implications of this gradient are discussed in depth in Section |HJ 



7.2. Age-Metallicity Relationships 



We will further investigate the nature of the changes in the spatial metallicity 
distribution with the help of the AMRs of each field. In spite of the large uncertainties 
of the individual age determinations presented in Section El the AMRs are useful to 
characterise the general trend of the run of metallicity with age. The AMR for each field 
is plotted in Figures |9] and [TOJ The many stars with the same age that appear both at the 
old and young limits result from the method degeneracy. As discussed in Paper II and in 



Carrera et al. 



(J2008bJ), our procedure to estimate ages saturates for values older than 10 



Gyr, so we can only give 10 Gyr as a lower limit for older stars. We assum e for them an age 



of 12.9 Gyr, the age of the oldest globular cluster observed in our galaxy ( ISalaris fc Weiss 



-18 - 



20021 . NGC 6426, ). In the case of the youngest stars in our sample, stellar evolution models 
indicate that the region in which we selected the observed stars is not populated by objects 
younger than 0.8 Gyr. Therefore, although equations used to estimate ages can formally 
lead to ages younger than this value, we assign them an age of 0.8 Gyr. 

All fields studied, independently of their distance to the center, show a rapid chemical 
enrichment reaching [Fe/H]>-1 in the first 2-3 Gyr. After this, the metallicity increased 
steadily but moderately until around 3 Gyr ago, thus increasing another ~0.5 dex in 
about 7 Gyr. This period corresponds wi th the age-gap observed in the clusters (e.g. 



Grocholski et al. 



2006 : 



Sharma et al. 



observed in the field (e.g. 



20 



Meschin et al 



Oj) an d with a period of lower star formation rate 



20 111 ). Finally, another episode of rapid chemical 



enrichment in the last ~3 Gyr is only observed in fields inner to ~7 kpc. In fact, we do not 
find stars younger than ~3 Gyr in our sample in the outermost fields despite the fact that 
these are among the most common on the innermost fields. The two pe riods of chemica 



enric hment at old and young ages are also observed in the clusters (e.g. 



Olszewski et al. 



199ll ). We refer the reader to Paper II for a detailed comparison among clusters, bar, and 
disk AMRs. In Paper II, we performed a x 2 test to check that the AMRs of the fields 
studied there are statistically indistinguishable among them. Due to the small number 
of stars confirmed as LMC member in the outermost fields studied here, a \ 2 test is not 
adequate in this case and we will restrict ourselves to a qualitative comparison among the 
AMRs of the different fields. To help with this comparison, we have plotted in Figures [9] 
andHnione of the chemical evolution models derived in Paper II that reproduces adequately 
the observed AMRs. This model was computed assuming a rate of gas outflow of 2.5 
(see Paper II for details j3- The AMRs of all fields are well reproduced by this model and 
therefore we feel confident that the chemical enrichment has proceded in a similar way in 



7 This model has been recomputed with the metallicity scale used in this paper. 
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all the studied fields. 

Inset panels of Figures 191 and [TU1 show the age distribution of the stars measured on each 
field, with and without taking into account the age uncertainty (solid line and histogram, 
respectively; see Paper II for details). In the innermost fields there are two peaks which 
correspond to the two periods of faster chemical enrichment. The amount of stars younger 
than 5 Gyr gradually decreases with increasing radius. Since the younger stars are also 
more metal-rich, the metallicity gradient observed in the outer disk may be explained by 
the gradual increase of the age of the youngest populations with radius. 



7.3. Kinematics of the LMC RGBs 



The kinematics of t 



sec 



van der Marel et ai- 



re LM C has been studied by many authors using different tracers 
20091 . for a review). There is agreement among workers that all 



tracers studied describe a similar cold disk kinematics but with different rotational velocities 



for each of them: 61 km s 1 for the carbon stars (Ivan der Marel et al. 



H I gas ( lOlsen fc Massey 



20071 ); and 107 km s 1 for red supergiants ( 



2002); 80 km s" 1 for 



Olsen fc Massey 



20071). 



It is also kno wn that the velocity dis persion, related to the disk scale height, increases from 



the youngest ( 



Olsen &; Massey 



( IBorissova et al. 



2006 



20071, red s upergiants ~9 km s 1 ,) to the oldest populations 



Minniti et al 



2003 



, RR Lyrae stars ~50 km s 1 ,). The fact that 



the oldest populations have the highest velocity dispersion suggests the possibility tha t 



they may be distribu ted in a halo similar to that of the MW (e.g. 



Minniti et al. 



Borissova et al. 



20061 ) 



20031 : 



Our sample of almost 500 RGB stars studied in an homogeneous way covers a 
large range of angular dist 
galactocentric radius (e.g. 



Zhao et al. 


2003; 


Cole et al. 


2005) 



20051 ). In Figure [TT] we have plotted 
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the velocities of the RGB stars studied here versus position angle for different galactocentric 
radius (filled circles). Our RGB sample is well re produced by the kinematical rotational 



disk model derived by 



van der Marel et al 



[TT1 the carbon stars data (dots) from 



(2002j) fo r carb on st ars (da s hed li ne). In Figure 



Kunkel et al. 



(119971 ) and lHardvi (12010h have also 



been plotted as comparison. The velocity dispersion of the RGB sample studied here 
(er=24.1±0.2 km s~ x ) is slightly larger than that of the carbon stars (q=20.2± 0.5 km s -1 ). 



Our r esult is in very good agreement with o ther analysis of LMC field RGB by 



(12003 



, a=24 km s 1 ,) and 



Cole et al 



(12005 



Zhao et al. 



a=24.7±0.4 km s -1 , 



Since RGBs have a wide range of metallicities and ages, we can compare the velocity 
dispersion of different ages/metallicities, which provides information about the scale height 
of different populations. In fact, if a classical halo formed by old, metal-poor, high-velocity 
dispersion objects existed in the LMC, a dependence between the velocity dispersion and 
metallicity might be observed. To investigate this we have compared the velocity dispersion 
in different metallicity bins. The mean velocity of each field was subtracted from the 
radial velocity of each star to eliminate the rotational component, as in Paper II. Then, 
we divided the metallicity range in 4 bins and computed the velocity dispersion of each 
of them by fitting a Gaussian to the velocity distribution. The number of objects and 
velocity dispersion in each metallicity bin are listed in Table It can be noted that the 
velocity dispersion increases as metallicity decreases, or age increases. This agrees with 
the fact found by other authors th at the velocity dispers i on increases from the you ngest 



population to the oldest ones (e.g. 



Olsen fc Massey 



2007 



van der Marel et al. 



20091). The 

velocity dispersion of our most metal-poor, and therefore older, bin is er=26.8±1.4 km s -1 
based on 45 objects. This value, si milar to that foun d in Paper II from a smaller sample, is 



lower than the dispersion found by 



Cole et al. 



(120051 ) for the most metal-poor RGBs in the 



bar, er=40.8±1.7 km s 1 based only on 18 stars. Although we believe that we are sampling 
the most metal-poor/oldest populations, we do not find evidences that they belong to a 
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hot halo. iBorissova et all (120061 ) (see also iMinniti et al.l 120031 ) found <r=50±2 km s -1 for a 
sample of RR Lyrae in the center of the LMC. These values are significantly higher than the 
value found here for the most metal-poor stars. These larger velocity dispersions have been 
quoted as evidence for the presence of a classical halo. However, this is still an open question 
since most RR Lyrae have an exponential surf ace density distribution with a similar scale 



lengt h as other well known tracers of the disk (lAlves 



2004 



Subramaniam fc Subramanian 



2009) , and only a small fraction, <1 0%, is likely to trace the extended halo of the LMC 



( ISubramaniam fc Subramanian 



2009|) 



The galactocentric distance covered by our sample also allows us to investigate the 
behaviour of the velocity dispersion as a function of radius. This is a key analysis to 
investigate the stability of the disk and, to our knowledge, this is the first time it is done 
in the LMC. As before, the velocity dispersion in four radial annuli has been computed by 
fitting a Gaussian to the rotational-free velocities of stars in each of them. The obtained 
values, together with the number of objects in each annulus, are listed in Table |6] and 
plotted as filled squares in Figure [12j The velocity dispersion decreases with radius with 
the exception of the outermost annulus studied. As comparison, we obtained in the same 



Cole et al. 



(120051 ) . The 



way the velocity dispersion of the Bar from the sample presented by 
Bar has a smaller dispersion than the innermost disk annulus analysed. Since our sample is 
not well distributed in angular distances, we have also calculated the velocity dispersion as 
a function of radius for the carbon stars sample described above (open circles). As in the 
case of the RGBs, the velocity dispersion decreases with radius if we exclude the central 
area. The behaviour in the outermost annulus studied may not be reliable because of the 
small number of objects studied at this galactocentric distance. The implications of these 
results are discussed in depth in the following section. 
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8. Discussion 

8.1. Age and Metallicity Gradients in the LMC 

The spatial distribution of our fields allows us to investigate the presence of age 
and/or metallicity gradients in the field populations of the LMC. In the 4 fields studied 
in Paper I, located at galactocentric distances from 3 to 7 kpc, we found that while the 
oldest populations are coeval in all fields, the age of the youngest component gradually 
increases with radius. With the new data, we confirm that the trend continues at larger 
galactocentric distances (~9 kpc). If we assume that the youngest stars were formed in-situ 
(e.g. there is no significant stellar migration), we are observing an outside- in quenching of 
star formation at recent times. 



This age gradient observed across the LMC disk is not accompanied by a gradient in 
the mean metallicity, which it is observed only for R>6 kpc. A gradient in the metallicity 
in the youngest objects cannot be excluded with our observations since they do not have 
counterparts in the red giant branch. However, the fact that the AMR is so similar at all 
radii points to a similar history of metal enrichment across the whole LMC disk. In this 
context, a gradient in the mean metallicity of different fields would be mainly related to a 
varying proportion of stars of different ages in dif ferent fields. A similar trend of metallicit y 



with radius has been found by 



Feast et al. 



(120 !0l ) for AGB stars (but see also 



Cioni 



2009|). 



Recent studie s in clusters a. 
(Figure 16 of ISharma et al- 



so sho w hints of a possible metallicity gradient in the outer disk 



2010l ). However, the small number of clusters at these outer 



distances prevents them from extracting further conclusions. 



A dynamical bar, as that obs erved in the 



of a disk within a certain radius (jFriedli et al. 



MC, might be able to mix the stellar content 



1994 



Zaritsky et al. 



19941 ) . Therefore this 



result might put a constrain to the mixing produced by the LMC bar. The age gradient 
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observed at R<6 kpc can be explained by the fact that the galaxy has not had enough time 
to mix populations younger than 1 Gyr. It is expected that the stellar content of the LMC 
is well mixed in a few crossing times. At a distance of 6 kpc, the crossing timej is ~2xl0 8 
yr which is in agreement with what observed. 

The mean metallicity of the most metal-rich objects in each field gradually decreases 
with the galactocentric distance from R~6 kpc. This trend is due to an increase of the 
age of the youngest population in each field with radius, because the chemical enrichment 



history has been the s ame over the 



in the M33 outer disk (IBarker et al 



jMC disk. A similar trend has been also observed 



20071 ) 



Taylor et al. 



(120051 ) studied the behaviour of 



color with radius in a sample of irregular and late-type spiral galaxies finding that these 
galaxies became redder with increasing radius which is also interpreted as the existence 
of an age gradient. This agrees with an outside-in disk evolution scenario proposed in 



Paper I for the 



(White k Frenk 



MC as opposed to t 



Trujillo fc Pohlen 



991 



Mo et al. 



that predicted by 



re inside-out scenario predicted by ACDM cosmology 



19981 ) and observed in more massive early-type spirals (e.g. 



200 5J). This o u tside- in disk evolution scenario is qualitatively similar to 



Stinson et al. 



(120091 ) f or dwarf galaxies an d observed in several Local 



Group dwarfs in the LCID 



the LMC, the 



p roject (e.g. 



Stinson et al. 



Hidalgo et al. 



20091 1 . However, for a galaxy like 



(120091 ) models predict an inside-out scenario. This may put 



constraints to the mechanisms that shape the radial distribution of stellar populations in 
galaxies of intermediate mass such as the LMC. 

It is well known that the LMC forms an interacting system with the SMC and the 
MW at the moment. However it is not clear if the sa me has happened in the past. New 



dynamical models, based on accurate proper motions (jKallivayalil et al. 



2006 



Piatek et al. 



^This value has been computed from the t c 



R/a (e.g. 



Binney fc Tremaine 



19871 ) 



using a velocity dispersion of <r=24.1 km s 1 derived in previous section. 
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20081 ). suggest that the Magella nic Clouds are un bound and they might be in their first 



passage around the Milky Way (IBesla et al. 



20071 ). Repeated encounters with one or more 



external systems might have erased, or at least weakened, any gradient. The results of this 
paper tend to support an evolution of the LMC in a more isolated environment. 



8.2. Kinematics of the RGB Field Populations 



Our RGB sample has rotational cold disk kinematics similar to that of carbon stars. 
The velocity dispersion is <r=24.1±0.2 km s _1 . It increases from most metal-rich/youngest 
stars, cr=20.9±1.2 km s -1 , to the mo st metal-poor /oldest ones, cr=26.8±1.4 km s _1 . Other 



tracers have a similar behaviour (see 



van der Marel et al. 



2009 



for a review). Therefore, 



the scale height of the LMC disk increases with time. However, even for the most 
metal-poor/oldest objects there are not evidences of the presence of a classical hot halo. In 
fact, recent works ind i cate that the LM C disk extends, at least, to a distance of 12-14 kpc 



(jMajewski et al. 



2009 



Saha et al. 



20101 ). 



The velocity dispersion of disk RGBs decreases with increasing galactocentric radius, 
except for those objects in the bar region, which have a lower velocity dispersion than those 
in the inner disk. Carbon stars have a similar trend. From pioneering works in the disk of 



spira 
(e.g. 



galaxies it is well known that the velocity dispersion decreases with increasing radius 



Bottema 



19931 ). The same trend has been also foun d in the 



galaxies studied un til now to our knowledge: NGC 2552 ( jSwaters 



only t wo disk irregular 



(IHunter et al. 



19991 ): and NGC 4449 



20051 ). This trend is easily explained by the decreases of stellar density with 



galactocentric dista nce if the scale height remain s constant with radius, as is observed in 



most disk galaxies ( Ivan der Kruit fc Searle 



19811 ). The fact that the velocity dispersion is 



lower in the bar than in the innermost disk i s a feature also observed in at least one half 
of barred spirals (e.g. IChung fe Bureau! 120041 ) . The velocity dispersion increase observed 
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in the outermost radius studied can be inter preted as a thic kening; in the external disk 
due to dynamical interactions with the MW (jWeinbergl |2000| ) . However the small numbe 
of objects confirmed as LMC members as this galactocentric distance prevents us from 
extracting further conclusions. 



9. Summary 

In this paper we have presented deep CMDs for 6 fields in the outer LMC disk, covering 
different distances, 5.9 kpc<R<9.2 kpc, and azimuths. Using CaT spectroscopy we have 
derived stellar metallicities for a sample of RGB stars in these fields, which have been 
analysed together with those stars studied in Paper II. In total, including those studied in 
Paper II, our sample has 497 disk RGB stars confirmed as LMC members from their radial 
velocity, covering a galactocentric radius from 3.0 to 9.5 kpc and different position angles. 
The main results of this paper are: 

• The comparison of the CMDs with isochrones shows that the age of the youngest 
populations gradually increases with radius. This is a continuation of the trend found 
in Paper I for inner fields. 

• The small number of stars confirmed as LMC members in the outermost fields (R>7 
kpc) does not allow a detailed analysis of the metallicity distributions of these fields. 
A statistical analysis shows that the 10th percentile of the metallicity distribution, 
which accounts for the most metal-poor populations, does not show a clear trend 
with radius within the errors. Therefore, we may conclude that the most metal-poor, 
and old, population is similar in all of them. On the contrary, the 90th percentile, 
indicative of the most metal-rich populations, remains constant in the inner 6 kpc, 
and gradually decreases from that point on. This result, which is independent of the 
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position angle, may be explained by the gradual increase of the age of the youngest 
populations with galactocentric distance. 

As in Paper II, we have estimated the ages of the observed stars from the combination 
of spectroscopic metallicities, and positions of the stars in the CMD. This allows us to 
obtain the AMR for each field. It is consistent with a chemical enrichment that has 
been similar in all studied fields as was found in Paper II for inner fields: i.e. coeval 
stars would have the same metallicity everywhere in the LMC disk. 

Our RGB sample has a rotational cold disk kinematics similar to that inferred 
from carbon stars, with a velocity dispersion of cr=24.1±0.2 km s _1 . The velocity 
dispersion increases from most metal-rich/youngest stars, <r=20.9±1.2 km s -1 , to 
the most metal-poor /ol dest ones, cr=26.8±l .4 km s _1 as has been observed for other 



kinematical tracers (see 



van der Marel et al. 



2009, for a review). However, even for 



the most metal-poor/oldest objects there is not evidence of the presence of a classical 
hot halo. 
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Fig. 1. — Spatial distribution of the fields studied in this paper (fille d squares) t ogeth er with 



those studied in Papers I and II (open circles). The field studied by lCole et al.l ( 120051 ) in the 



bar is also marked (filled triang le). As reference, the isopleths of the RGB and AGB stars 



derived by 



van der Marell (120011 ) have been superimposed with semimajor-axis of 1°, 2°, 3 C 



4°, 6°, 8°, and 10°. The position of the SMC is also shown in the bottom-right corner. Right 
ascension and declination are labeled on top and right axis, respectively. 
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Fig. 2. — {(B — I) ,Mj] CMDs for the eastern (top) and western (bottom) fields. Isochrones 



with a ges and metallicities as labeled, and a zero-age horizontal branch by iPietrinferni et al. 



( 120041 ) have been superimposed. The identification of each field together with its galacto- 



centric distance, in kpc, are labeled in top-left corner. 
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Fig. 3. — [(V — I)o,Mj] CMDs for the northern fields. As in Figured isochrones have been 
superimposed. The identification of each field together with its galactocentric distance, in 
kpc, are labeled in top-left corner. 
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Fig. 4. — Top panels: [Fe/H] vs. W'y (left) and W'j (right). The solid lines are the best 
linear fits to the data. Dashed lines represent the 90% confidence level band of the fits. Open 
circles are clusters younger than 4 Gyr. Bottom panels: the residuals of the linear fits. Note 
that the W errors are smaller that the size of points in several cases. 



-37- 



= 0.2 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
7 Paper II ~ 


i i i i i i i i i i i i i i i i i 

: S09 •-. 


i , , , , i , , , , i , , , , i , 


. . ..rf .1' . . 
• • ■ • • 

i , , , " i , ;", ; i , t , , i , 



-3 -2 -1 -3 -2 -1 

[Fe/H] Here 



Fig. 5. — Differences between metalli cities derived here with th ose obtained in Paper II (left 
panel) and those calculated following IStarkenburg et al.l (1201(1 right panel) as a function of 
the values used in this paper. For this comparison we used only stars in the northern fields. 
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Fig. 6. — Metallicity distributions for the northern fields. The innermost field is on the 
topdeft and the outermost one on the bottom-right. The galactocentric distances in kpc of 
each field are labeled on the top-left corner. The four innermost fields were the subject of 
Paper II, for which metallicities have been recalculated in the new scale defined in Section [5j 
The errorbars are the square root of the number of stars in each bin. 
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Fig. 7. — The same as Figure El for the eastern (top) and western (bottom) fields. 
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Fig. 8. — Metallicity at the 10th (metal-poor; green filled symbols) and 90th (metal-rich; 
red empty symbols) percentile and mean (black starred symbols) values of the metallicity 
distribution of each field as a function of projected radius. See text for details. 
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Fig. 9. — Age-metallicity relationships for the northern fields in our sample. One of the 
chemical enrichment models derived in Paper II has been overplotted for comparison (see 
text for details). The galactocentric distances of each field are labeled in top-right corner, in 
kpc. Inset panels show the age distribution computed by taking (solid line) and by not taking 
(histogram) into account the age determination uncertainties (see Paper II for details). Age 
error are given in the top strip while metallicity errors are shown in the right one. 
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Fig. 10. — The same as Figure El for eastern (top) and western (bottom) fields in our sample. 
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Fig. 11. — Mean and dispersion of the velocity distributions (large dots with errorbars) as 
a function of the position angle (note that each distribution is plotted twice) at different 



angular distances (from top-left to bottom-r ight). As comparison, car 



Don stars data have 



Kunkel et al 
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are the predictions of the kinematical model computed by 
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Fig. 12. — Radial velocity dispersion as a function of radius for the RGBs (filled squares) 
and for the carbon stars (open circles). With the exception of the bar region, both show a 
decrease of the velocity dispersion with radius up to a distance of 6 kpc. From there off, an 
increase of velocity dispersion is observed although the small number of objects studied at 
these galactocentric distances do not allow us a reliable discussion. 
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Table 1. Observed Fields 



ID 




«2000 


^2000 


p a 


P. A. a 


E(B-V) 


p b 


R(kpc) c 


LMC0512- 


•6648 


05 h 12 m 


-66°48' 


3?4 


333?0 


0.10 


3? 5 


3.0 


LMC0514- 


•6503 


05 h U m 


-65°03' 


5?0 


343?3 


0.05 


5?0 


4.4 


LMC0513- 


•6333 




-63°33' 


6?5 


345?4 


0.04 


6?5 


5.7 


LMC0422- 


•7100 


04 /l 22 m 


-71°00' 


5?6 


250?7 


0.08 


6?8 


5.9 


LMC0513- 


•6159 




-61°59' 


8?0 


347?6 


0.03 


8?0 


7.1 


LMC0649- 


■7110 


06 /l 49 m 


-71°10' 


6?9 


110?3 


0.13 


8?0 


7.1 


LMC0512- 


•6100 




-61°00' 


9?0 


347?8 


0.03 


9?0 


7.9 


LMC0356- 


■7100 


03 h 56 m 


-71°00' 


7? 7 


250?9 


0.09 


9?3 


8.3 


LMC0511- 


•6000 




-60°00' 


10?0 


348?0 


0.03 


10?0 


8.8 


LMC0714- 


■7100 


07 h U m 


-71°00' 


8?9 


109?7 


0.19 


10?4 


9.2 



a Angular coordinates over the celestial sphere (see text for details). 
b De-projected radius calculated assuming i=3A°.7 and 6=189?3. 
b Distance in kpc calculated as R[kpc] = D x tan(d[ ]) assuming D=51 kpc. 



-46- 



Table 2. LMC RGB Stars Observed. Only a portion of this table is shown here to 
demonstrate its form and content. A machine- readable version of the full table is available. 



«2ooo ^2000 SCa UaCa B or V I V r (km s x ) a Vr (km s r ) Comments 



04:23:02.04 


-70:58:51.0 


8.218 


0.262 


19.9325 


17.9836 


239.6648 


2.880 


04:22:52.95 


-70:57:00.3 


8.605 


0.147 


19.4462 


16.8805 


251.9283 


1.260 


04:23:25.27 


-70:56:16.2 


7.839 


0.120 


18.9702 


16.4791 


241.2483 


1.494 


04:22:21.24 


-70:55:42.2 


7.565 


0.203 


19.8558 


17.7160 


225.6490 


2.196 


04:23:00.06 


-70:54:18.2 


8.422 


0.102 


19.3202 


16.9670 


214.9008 


1.587 



Low S/N 
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Table 3. Cluster Sample 



Cluster WV W'j [Fe/H] CG97 [Fe/H] C09 Reference 



NGC 104 


6 


.94±0. 


,01 


6. 


,23±0 


.01 


-0 


.67±0. 


.05 


-0.76±0.02 


1 


NGC 188 


8. 


17±0. 


.07 


7. 


.27±0. 


.08 


-0. 


07±0. 


.04 


+0.01±0.08 


2 


NGC 288 


5. 


.51±0. 


.01 


5 


,04±0. 


.01 


-1. 


,07±0 


.03 


-1.32±0.02 


1 


NGC 362 


6. 


,01±0. 


.01 








-1. 


09±0. 


.03 


-1.30±0.04 


1 


NGC 1851 


5. 


,94±0. 


03 


5 


,24±0. 


.04 


-1. 


03±0 


.06 


-1.18±0.08 


1 


NGC 1904 


4 


,91±0. 


.03 








-1. 


37±0. 


.05 


-1.58±0.02 


1 


Be 20 


6 


.86±0. 


03 


6. 


.39±0. 


03 


-0. 


49±0. 


.05 


-0.49±0.05 


3 


NGC 2141 


8. 


,33±0. 


.01 


7. 


,67±0. 


02 


-0. 


14±0. 


.05 


+0.00±0.16 


4 


Cr 110 


8. 


,21±0. 


.04 


7. 


,74±0. 


06 








+0.03±0.02 


5 


NGC 2298 


3 


75±0. 


03 


3 


,09±0. 


03 


-1. 


74±0. 


.04 


-1.96±0.04 


1 


Melote 66 


7 


,69±0. 


03 


6. 


,90±0, 


03 


-0. 


38±0. 


.06 


-0.33±0.03 


3 


Be 39 


8. 


.21±0. 


.04 


7. 


.27±0. 


06 








-0.21±0.01 


6 


NGC 2682 


8. 


,24±0. 


01 


7. 


,48±0. 


01 


-0. 


03±0 


.03 


+0.05±0.02 


5 


NGC 3201 


5. 


.46±0. 


03 


4 


.76±0. 


02 


-1. 


23±0 


.05 


-1.51±0.02 


1 


NGC 4590 


2. 


,84±0. 


.02 


2 


,19±0. 


06 


-1. 


,99±0 


.03 


-2.27±0.04 


1 


NGC 5927 


7. 


,81±0. 


.01 


6 


,92±0. 


01 








-0.29±0.07 


1 


NGC 6352 


7 


10±0. 


.01 


6 


,31±0. 


01 


-0. 


64±0 


.06 


-0.62±0.05 


1 


NGC 6528 


8. 


.28±0. 


.04 


7. 


.58±0. 


.04 


-0. 


17±0. 


.02 


+0.07±0.08 


1 


NGC 6681 


5. 


.05±0. 


03 


4. 


,49±0. 


.07 


-1. 


35±0. 


.03 


-1.62±0.08 


1 


NGC 6705 


8. 


.95±0. 


.07 


8. 


,28±0. 


.12 


+0.07±0.05 


+0.07±0.07 


7 


NGC 6715 


4. 


.86±0. 


03 


4 


,30±0. 


03 


-1.25±0.06 


-1.44±0.07 


1 
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Table 3 — Continued 



Cluster 


wv 




[Fe/H] CG97 


[Fe/H] C09 


Reference 


NGC 6791 


9.78±0.09 


8.77±0.09 


+0.47±0.04 


+0.47±0.04 


8 


NGC 6819 


8.41±0.04 


7.64±0.04 


+0.07±0.03 


+0.09±0.03 


9 


NGC 7078 


2.78±0.01 


2.18±0.01 


-2.12±0.04 


-2.33±0.02 


1 


NGC 7789 


8.47±0.02 


7.61±0.02 


-0.04±0.05 


-0.04±0.05 


10 



References. 


- (1) Carretta et al. ( 


200 


3); 


'2) 


Randich et a 


1. (200? 


); (3) 


Yone et al. 


C 


2005): (4) J 


acobson et al 


• 


2009 


); 


'5) Pancino et al. 


< 


2010): 


(6) Friel et al 


(2010); (7) 


Gonzalez & Wallerstein 


(2000 N 


(8) ( 


Carretta et al. 


(2007); (9) 


Braeaelia et al. 


(2001); (10) 


Tautvaisiene et al 


. (2005 


)• 



Table 4. Mean, standard deviation, 90th and 10th percentile values of the metallicity 

distribution of each field. 



ID R(kpc) ([Fe/H]) a [Fe/H] [Fe/R] 90th [Fe/R] 10th N* 



Bar 




0.3 


-0.40±0.07 


n or i n 

0.35±0. 


07 


-0. 


25±0. 


.12 


-1 


,25±0. 


on 

.30 


373 


LMC0512- 


•6648 


3.0 


-0.58±0.09 


0.47±0. 


07 


-0. 


.09±0. 


.09 


-1 


.38±0. 


.28 


130 


LMC0514- 


■6503 


4.4 


-0.64±0.10 


0.54±0^ 


11 


-0. 


.09±0. 


10 


-1 


,71±0. 


.44 


110 


LMC0513- 


■6333 


5.7 


-0.54±0.11 


0.48±0. 


.07 


-0. 


,01±0, 


.08 


-1 


,57±0. 


.17 


69 


LMC0422- 


■7100 


5.9 


-0.43±0.09 


0.39±0. 


.14 


-0. 


,09±0. 


.07 


-1 


,09±0. 


.48 


63 


LMC0513- 


•6159 


7.1 


-1.03±0.13 


0.64±0. 


.10 


-0. 


,27±0. 


.16 


-2 


17±0. 


.25 


65 


LMC0649- 


■7110 


7.1 


-0.74±0.14 


0.62±0. 


.15 


-0. 


,10±0. 


.15 


-1 


,79±0. 


.14 


24 


LMC0512- 


•6100 


7.9 


-1.04±0.30 


0.70±0^ 


.25 


-0. 


.42±0. 


.15 


-1 


.96±0. 


.47 


12 


LMC0356- 


■7100 


8.3 


-0.87±0.23 


0.56±0. 


.16 


-0. 


.48±0. 


.08 


-2. 


,07±0. 


.25 


12 


LMC0511- 


-6000 


8.8 


-0.99±0.37 


0.20±0. 


.10 


-0. 


.71±0. 


.10 


-1 


.63±0 


.47 


6 


LMC0714- 


■7100 


9.2 


-0.85±0.22 


0.35±0. 


23 


-0. 


,71±0. 


.07 


-1 


.75±0 


.14 


7 



Table 5. Velocity dispersion for each metallicity bin. 



[Fe/H] 


Age (Gyr) 


N* 


<jy(kms v ) 


-0.25 to +0.5 


<2 


117 


20.9±1.2 


-0.75 to -0.25 


2-6 


212 


21.8±0.9 


-1.5 to -0.75 


6-10 


122 


26.2±1.1 


-3.5 to -1.5 


>10 


45 


26.8±1.4 
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Table 6. Velocity dispersion clS db function of radius. 



R (kpc) 


N* 


av(kms 1 ) 


Bar 


374 


25.4±0.2 a 


2.5-3.5 


130 


28.6±1.5 


3.5-5 


172 


23.6±0.4 


5-7 


138 


21.1±1.2 


7-10° 


56 


24.0±1.1 



a Valuc obtained from Cole et al 



i 20051 ) data following the same pro- 
cedure than for our sample. This 
value is similar to that obtained in 
the original work of 24.7±0.4 km 
=-1 



